Stable Patterns of Membrane Domains at Corrugated Substrates 



OO 
O 

o 

(N 



O 



I 

O 

o 



> 

m 
i> 
cn 

o 

00 

o 



X 



Bartosz Rozycki, Thomas R. Weikl and Reinhard Lipowsky 
Max Planck Institute of Colloids and Interfaces, Science Park Golm, 144^4 Potsdam, Germany 

(Dated: February 4, 2008) 

Multi-component membranes such as ternary mixtures of lipids and cholesterol can exhibit coex- 
istence regions between two liquid phases. When such membranes adhere to a corrugated substrate, 
the phase separation process strongly depends on the interplay between substrate topography, bend- 
ing rigidities, and line tension of the membrane domains as we show theoretically via energy mini- 
mization and Monte Carlo simulations. For sufficiently large bending rigidity contrast between the 
two membrane phases, the corrugated substrate truncates the phase separation process and leads 
to a stable pattern of membrane domains. Our theory is consistent with recent experimental obser- 
vations and provides a possible control mechanism for domain patterns in biological membranes. 

PACS numbers: 87.16.D-, 87.16. dt, 64.75. St 



Introduction - Biomimetic and biological membranes 
are 2-dimensional liquids, in which the lipid molecules 
undergo fast lateral diffusion. In general, the interac- 
tions between the different lipid species may lead to the 
formation of intramembrane domains with distinct lipid 
composition. For giant vesicles prepared from three- 
component membranes, such domains can be directly 
observed by optical microscopy [l|, 0, H, 0j Si • These ex- 
periments have confirmed the process of domain-induced 
budding as predicted theoretically @. The latter shape 
transformation provides direct evidence for the line ten- 
sion of the domain boundaries, which was found to vary 
between 10~^^ and fO~^'* N for different compositions 
and temperatures 0] and must vanish at the critical 
demixing or consolute point of the membrane mixture 

S- 

The three-component membranes studied in [l|, [3, 0, 0j 
consisted of a saturated lipid such as sphingomyelin, 
an unsaturated phospholipid, and cholesterol. The same 
lipid mixtures were proposed to form domains in cell 
membranes 0, @] but direct imaging of these latter do- 
mains or 'rafts' has turned out to be difficult. In con- 
trast, it is now well established that the three-component 
model membranes exhibit two-phase coexistence regions 
in which the membrane components separate into liquid- 
ordered (Co) and liquid-disordered {Cd) domains that are 
enriched in the saturated and unsaturated lipids, respec- 
tively. 

The micrometer-size domains observed in three- 
component membranes must arise from the growth and 
coarsening of much smaller domains, a process which acts 
to reduce the line energy of the domain boundaries be- 
tween the Co and Cd domains. This coarsening process 
can be modified by constraints on the membrane shape 
as observed in recent experiments [^, , in which the 
membranes adhered to a substrate surface. For a planar 
substrate, the membrane domains should again undergo 
the same coarsening process as for giant vesicles. For a 
corrugated substrate surface as in Fig. [U on the other 
hand, the domains tend to form certain domain patterns 



as observed in [^, [lo| . One important and open ques- 
tion is whether these domain patterns are metastable or 
represent the true equilibrium states of the membranes. 

In this letter, we will address and answer this question. 
First, we will clarify the basic mechanism underlying the 
observed domain patterns and show that these patterns 
arise from the competition between the line tension. A, of 
the domain boundaries and the bending rigidities Kq and 
Kd of the ordered and disordered liquid phases, where 
we take Kq > as in The bending rigidity con- 
trast An = Ko — Kd and the line tension A define the 
elastic length scale ^ci = Ak/A which must be compared 
with another length scale, ^to, that is determined by the 
surface topography alone. If the elastic length scale ^ei 
exceeds the topographical length scale ^to, the pattern 
with many striped domains is predicted to be globally 
stable and to represent the true equilibrium state of the 




FIG. 1: Multi-component membrane adhering to a corrugated 
substrate surface - In this example, the topography is peri- 
odic in one direction and characterized by the curvature radii 
±i? of the cylindrical ridges and valleys, the width L of the 
flat surface segments, and the tilt angle (j). The membrane 
contains liquid ordered {£o, gray) and liquid disordered {£d, 
white) domains. The £o phase, which has a higher bending 
rigidity than the jCd phase, tends to avoid the curved mem- 
brane regions along the ridges and valleys. 
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membrane. For ^oi < Cto, on the other hand, the equi- 
hbrium state consists of one large Co and one large Ld 
domain, both of which cover many ridges and valleys of 
the surface topography. 
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Our letter is organized as follows. We first define our 
theoretical model in terms of the different energies of a 
multi-component membrane that forms two types of in- 
tramembrane domains and adheres to a topographically 
structured surface. Using this model, we then determine 
the phase diagram both for the grand-canonical ensem- 
ble, in which the membrane is in contact with reservoirs 
for the different lipid molecules, and for the canonical en- 
semble, in which the membrane has a constant number of 
lipid molecules. Both phase diagrams are quantitatively 
confirmed by extensive Monte Carlo simulations, which 
provide additional insight into the anisotropic coarsening 
process of the membrane domains. 

Theoretical Model - The system under consideration 
is a multicomponent lipid membrane, which strongly ad- 
heres to a topographically structured substrate as in 
Fig. [TJ If the membrane is quenched into (or pre- 
pared within) the Co^Cd coexistence region, it forms two 
types of membrane domains corresponding to these two 
phases and characterized by bending rigidities Kq and Kd- 
The bending energy of the membrane is then given by 
= J^^ dA 2koM2 + J^^ dA 2KdM^ where M denotes 
the mean curvature of the adhering membrane and the 
two surface integrals extend over the total area of the Co 
and £-d domains, respectively. 

The energy of the domain boundary separating the two 
lipid phases in the membrane is Su — Jg^ dl A, which 
depends on the line tension A > and the total length 
Jg^ dl of the domain boundaries. The competition of 
bending energy £bc and line tension energy Sn was first 
studied in [6| for multi-domain vesicles. For adhering 
membranes as considered here, we must also include the 
adhesion energy of the membrane which has the form 
= -J^^dAWo - li^^dAWd, where Wo > and 
Wd > denote the adhesion energy per unit area of the 
Co and Cd membrane segments, respectively. Combin- 
ing all three energy contributions, the total energy of the 
membrane conformation can be written as 



dA (2 Ak - AW) 



dlX 



(1) 



with the bending rigidity contrast Ak, = Kq — Hd and the 
adhesion energy contrast AW = Wo — Wd- Note that 
AW > corresponds to a substrate that attracts the 
ordered Co phase more strongly than the disordered Cd 
phase. 

In the following, we will focus on the surface topogra- 
phy displayed in Fig. [TJ In this case, the expression ([!]) 
for the total energy simplifies since the mean curvature 
M assumes only three distinct values: M = l/(2i?) along 
the cylindrical ridges, M = — l/(2i?) along the cylindri- 
cal valleys, and M = along the planar surface segments. 



Phase diagram: Grand-canonical ensemble - From the 
theoretical point of view, it is convenient to first consider 
the grand-canonical ensemble, in which the membrane 
is coupled to reservoirs for the different lipid molecules 
and the area fractions of the two types of membrane do- 
mains can adjust freely. In this case, minimization of 
the total energy £ as given by H]) leads to three differ- 
ent membrane phases. In phases (I) and (II), the whole 
membrane segment that adheres to the substrate is in 
the Cd phase and in the Co phase, respectively. In phase 
(III), on the other hand, the two phases separate in such 
a way that the membrane stripes that adhere to the flat 
pieces of the substrate are in the more rigid Co phase 
whereas the membrane parts that adhere to the ridges 
and valleys of the substrate stay in the Cd phase. More 
precisely, the membrane attains phase (I) if 
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Inspection of the inequalities ©-([l]) shows that the 
overall phase diagram is determined by three dimen- 
sionless parameters as given by (jiAn/XR, LAW/X, and 
2(j)R/L. The first two parameters are proportional to the 
bending rigidity contrast An and to the adhesion energy 
contrast AW, respectively, and inversely proportional to 
the line tension A. The third parameter 2(j)R/L depends 
only on the topography of the substrate surface, com- 
pare Fig. [TJ The phase diagram corresponding to fixed 
topography with 24>R/L = 1/3 is displayed in Fig. [5J 

The three regimes in Fig. [2l may be distinguished by 
the area fraction Xo of the Co domains. This fraction 
is defined via Xo = Ao/{Ao + Ad), where Ao and Ad 
represent the total areas of the Co and Cd domains. In 
regimes (I), (II), and (III), one has Xo = 0, Xo — I, 
and Xo = 1/(1 + 2(t)R/L), respectively, which implies 
first order transitions between these regimes. The three 
phase boundaries, which correspond to the dashed lines 
m Fig. [2 meet in a triple point. This point is located at 
2A = L AW and 2A (1 + 2(/)i?/L) = Ak/R. 

Phase diagram: Canonical ensemble - Next, let us con- 
sider an adhering membrane with a certain, fixed lipid 
composition. In equilibrium, the area fraction Xo of the 
Co domains then attains a fixed value as well, and the 
adhesion energy contrast AW now plays the role of a 
Lagrange multiplier in ([TJ). The corresponding phase di- 
agram is shown in Fig. [3J as a function of area fraction 
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FIG. 2: Phase diagram for the grand-canonical ensemble as 
a function of the two dimensionless parameters (j>AK/\R and 
I/AW^/A for 2(j}R/L = 1/3. The dashed phase boundaries have 
been obtained by minimization of the total membrane energy 
£ as given by ((1} , the data points by Monte Carlo simulations. 
Regimes (I) and (II) correspond to uniform membranes in the 
Cd and Co phase, respectively. Regime (III) is determined by 
the inequalities ([4} and corresponds to a striped membrane in 
which the membrane domains that adhere to the flat pieces of 
the substrate are in the more rigid Co phase while the domains 
that adhere to the ridges and valleys stay in the Cd phase. 



Xo and dimensionless rigidity contrast (j)AK/XR. The 
three regimes (I) , (II) , and (III) of the grand-canonical 
phase diagram in Fig. [5] are mapped onto three vertical 
lines with constant Xo in the canonical phase diagram, 
see Fig. |3l The coexistence lines in the grand-canonical 
phase diagram, on the other hand, are now mapped into 
three coexistence regions of the canonical phase diagram 
as shown in Fig. [3l 
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FIG. 3: Phase diagram in the canonical ensemble as a func- 
tion of area fraction Xo and dimensionless rigidity contrast 
(fiAn/XR. The three phases (I), (II), and (III) are now rep- 
resented by three vertical lines with Xo = 0, Xo — 1 and 
Xo ~ L/{2(f)R + L), respectively. For intermediate values of 
the area fraction Xo, one has coexistence regions in which 
two of the three phases coexist. The horizontal dotted line 
with (pAn/XR — 2{2<j)R + L)/ L represents transitions from a 
membrane with many striped domains to a membrane with 
two large domains. 

The horizontal dotted line in Fig. [3] represents transi- 



tions from a membrane with many striped domains to a 
membrane with two large domains. This line is given by 
4>Ak/XR = 4>^oi/R = 2{2(j)R + L)/L. The latter relation 
implies that the elastic length scale ^oi = Ak/A must be 
compared with the topographical length scale 



ito = 2R{2ct>R + L)/{4>L)., 



(5) 



which depends only on the substrate topography, see 
Fig. [TJ For ^ei > ^to, the bending rigidity contrast Ak 
dominates over the line tension A and the membrane con- 
tains many striped Co and Cd domains as confirmed by 
Monte Carlo simulations, see Fig.m For ^ci < Cto? on the 
other hand, corresponding to large line tension and/or 
small bending rigidity contrast, the membrane will com- 
pletely phase separate into two large domains. 

Monte Carlo simulations - The phase diagram ob- 
tained by energy minimization has been confirmed by 
Monte Carlo (MC) simulations. To perform such simula- 
tions, the membrane surface is divided up into square 
patches of side length a. This leads to a square lat- 
tice that is composed of parallel stripes of length L||. 
These stripes have alternating width L and 2(f)R, com- 
pare Fig.[l] where the stripes with width 2(/)i? correspond 
to the curved ridges and valleys of the substrate surface 
and will be denoted by C. Each patch is labeled by a 
pair of integer numbers i — {ix,iy) with 1 < ix ^ 
and 1 < iy < Ny. Periodic boundary conditions are im- 
posed in both directions. 

On each patch i, we place the occupation number 
with Tii = Q and = 1 if the square patch i is occupied 
by an Cd and an Co domain, respectively. The discretized 
membrane energy Sdm is then given by 



^dis 
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corresponding to the bending, adhesion, and line energy 
as in ([1]). The first term contains a summation over all 
patches i that are contained in the curvature stripes C 
while the last term contains a summation over all pairs of 
nearest neighbor patches {ij). The coupling constant Aq 
is related to the line tension A via \a/T = Xoa/T — 
ln[(l -I- e~^''°/^)/(l - e-^f"/^)] with temperature T in 
energy units. 

For the grand-canonical ensemble, we determined the 
equilibrium states of the discretized model as given by ([6]) 
using Glauber dynamics . During each move of this 
MC algorithm, a membrane patch i is chosen randomly 
and the value of corresponding variable Ui is changed 
from to 1 or from 1 to 0. This trial move is then ac- 
cepted according to the Metropolis criterion. In this way, 
we studied the stability of different domain patterns and 
determined the critical nucleation size of the domains 
as a function of the model parameters. Discontinuous 
changes in the critical nucleation size of the domains in- 
dicate transitions between the homogeneous membrane 
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phases (I) or (II) and the striped membrane phase (III). 
Therefore, the domain stabihty analysis based on the MC 
simulations allows us to determine the phase boundaries 
where parallel stripes of the Lo and Cii phase coexist in 
the membrane, see Fig. [5J For small lattice parameter a 
with a <C min{2^i?, iy}, our simulation results are in a 
good agreement with the mean field predictions ([4]) . The 
MC data included in Fig. [2] were obtained for i/a = 120. 




FIG. 4: Two different snapshots of membrane domains as ob- 
tained from canonical Monte Carlo simulations for area frac- 
tion Xo = 0.4 after (top) 10^ and (bottom) 10^ MC moves 
per patch. The black and white regions correspond to Co and 
Ld patches, respectively, [l^ 

We also studied the process of domain formation 
within the canonical ensemble, in which the number 
M = '^i of membrane patches in the Co phase is 
kept constant. This constraint was incorporated via dif- 
fusive Kawasaki dynamics fl3| , in which two neighboring 
patches i and j with = 1 and Uj = are chosen ran- 
domly and the values of ni and Uj are swapped accord- 
ing to the Metropolis criterion. Starting from a random 
initial configuration {nj}, we then observed anisotropic 
coarsening processes as illustrated in Fig. [H In this ex- 
ample, the model parameters fulfill the inequality (O cor- 
responding to coexistence between phases (I) and (III). 
As the system evolves in time, the more rigid patches 
start to aggregate in the flat membrane stripes and to 
form small domains of irregular shapes. Some of these 
domains shrink and, eventually, vanish whereas other do- 
mains grow and become elongated parallel to the flat 
membrane stripes. Since the diffusion-limited exchange 
process across the ridges and valleys between the flat 
stripes, is rather slow, we performed additional simula- 



tions starting from the state of complete phase separation 
and confirmed that this state is stable apart from shape 
fluctuations of the 1-dimensional interface between the 
coexisting phases (I) and (III). 

Summary - We have theoretically studied domain pat- 
terns in multi-component membranes supported on cor- 
rugated substrates as in Fig. [TJ We showed that the pro- 
cess of pattern formation is governed by the competition 
between bending rigidity contrast and line tension. The 
corresponding phase diagrams are displayed in Fig. [2] and 
in Fig. [3] for the grand-canonical and canonical ensem- 
ble, respectively. The stability of substrate-induced do- 
main patterns can be understood in terms of two length 
scales, the elastic length ^oi = An/X and the topograph- 
ical length ^to as given by (O. The substrate induced 
patterns with many striped membrane domains repre- 
sent the true equilibrium states of the system provided 
Cci > £,to- For the coexistence regions of the Co and 
Cd phases as studied experimentally, the elastic length 
scale ^ei is of the order of 100 nm as follows from the 
experimental data in [2] . This implies that both regimes 
^ci > Cto and ^oi < Cto should be experimentally accessi- 
ble. 
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